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Abstract 

The paper gives a quick account of the simplest cases of the Hitchin integrable sys- 
tems and of the Knizhnik-Zamolodchikov-Bernard connection at genus 0, 1 and 2. In 
particular, we construct the action- angle variables of the genus 2 Hitchin system with 
group SL 2 by exploiting its relation to the classical Neumann integrable systems. 



1 Hitchin systems 



As was realized by Hitchin in [16|, a large family of integrable systems may be obtained 
by a symplectic reduction of a chiral 2-dimensional gauge theory. Let E denote a closed 
Riemann surface of genus g and let G be a complex Lie group which we shall assume 
simple, connected and simply connected. We shall denote by A the space of Lie(G ! )-valued 
0,1-gauge fields^ A = A^dz on E. Hitchin's construction associates to E and G an 
integrable system obtained by a symplectic reduction of the infinite-dimensional complex 
symplectic manifold T*A of pairs (A,$>) where = <3? 2 dz is a Lie(G)-valued 1,0-Higgs 
field. The holomorphic symplectic form on T*A is 

tr 5$ 5A (1.1) 



where tr stands for the Killing form on Lie(G) normalized so that tr<f> 2 = 2 for the long 
roots 4>. The local gauge transformations h £ Q = Map(Ti, G) act on T*A by 

A i — > h A = hAhT 1 + hdh~ 1 , $i — > h $ = h^h' 1 (1.2) 

preserving the symplectic form. The corresponding moment map fi : T*A — > Lie(Q)* = 
A 2 (E) (g> Lie(G) takes the form 

H(A,$) = d$ + A<5> + $A (1.3) 

The symplectic reduction gives the reduced phase space 

V = M _1 ({0})/5 (1-4) 

1 one may work in a fixed smoothness class and use the Sobolev norms to define topology in A 



with the symplectic structure induced from that of T*A. V may be identified with the 
complex cotangent bundle T*J\f to the orbit space M = AjQ and M, in turn, with the 
moduli space of holomorphic G-bundles on E. More precisely, care should be taken to avoid 
non-generic bad orbits in order to obtain tractable orbit spaces. This may be done by con- 
sidering only gauge fields A leading to stable G-bundles forming smooth moduli space M s 
or those leading to semi-stable bundles giving rise to a, generally singular, compactification 
M ss of M s - In what follows we shall be somewhat cavalier about such details. 

The Hitchin system has V as its phase space. Its Hamiltonians are obtained the fol- 
lowing way. Let p be a homogeneous Ad-invariant polynomial on Lie(G) of degree d p . 
Then 

hp(A,&) = p(^=p(^ z )(dz) d " (1.5) 

defines a dp-differential on S which is holomorphic if /j,(A, <£) = 0. Since h p is constant on 
the orbits of Q, it descends to the reduced phase space: 

h p : V — ► H°(K d »). (1.6) 

Here K stands for the canonical bundle (of covectors oc dz) and H°(K dp ) is the (finite- 
dimensional) vector space of the holomorphic cip-differentials on E. The (components 
of) hp Poisson-commute (they Poisson-commute already as functions on T*A since they 
depend only on the "momenta" <£). The point of Hitchin's construction is that, by taking 
a complete system of polynomials p, one obtains on V a complete system of Hamiltonians 
in involution. For the matrix groups, the values of Hamiltoniens h p at a point of V may 
be encoded in the spectral curve C obtained by solving the characteristic equation 

det($-f) = (1.7) 

for £ E K. The spectral curve of the eigenvalues £ is a ramified cover of E. The correspond- 
ing eigenspaces of <3? form then a holomorphic line bundle over C belonging to a subspace 
of the Jacobian of C on which the Hamiltonians h p induce linear flows. 

For the quadratic polynomial p2 = |tr, the map h P2 takes values in the space of holo- 
morphic quadratic differentials H°(K 2 ). This is the space cotangent to the moduli space 
of complex curves E. Variations of the complex structure of E are described by Beltrami 
differentials 5fi = b[i% d z dz such that z' = z + Sz with dzSz = gives new complex 
coordinates. The Beltrami differentials b\i may be paired with quadratic differentials j3 by 

{PM ^ I f35fi. (1.8) 

The differentials Sfi = d(5£), where <5£ is a vector field on E, describe variations of the 
complex structure due to diffeomorphisms of E and they pair to zero with (3. The quotient 
space H l (K~ l ) of differentials 5fi modulo 8(5£) is the tangent space to the moduli space of 
curves E and H°(K 2 ) is its dual. The pairing (|1.8| ) defines then for each [5fj] G i7 1 (i^ _1 ) 
a Hamiltonian 

hsfi = / hp 2 5/i. (1.9) 



2 



The Hamiltonians hgn Poisson-commute for different Sfi. 



Hitchin's construction possesses a natural generalization [pO]1 [j5|| Iff] . Let x n £ E be 
a finite family of distinct points in E and O n a family of (co) adjoint orbits in Lie(G)* = 
Lie(G). 

O = {J2 x nS Xn | \ x eO n }, (1.10) 

n 

where 5 X stands for the Dirac delta measure at x, forms a coadjoint orbit of the group Q of 
local gauge transformations. In the symplectic reduction we may replace definition ( jL4| ) 
with 

V a = ir\{0})/G - fi'HYl X n$x n )/G\,x (i.n) 

where G\, x is the subgroup of Q fixing J2 Ki$x„ ■ Upon restriction to properly defined stable 



pairs (A, <&), V gives a smooth space with a semi-stable compactification [20]. Its second 



representation in ( 1.11| ) allows to equip V with the symplectic structure inherited from 



T*A. The reduced Hamiltonians h p take now values in H°(K dp (d p J2 x n)), i-e. define 
meromorphic dp-differentials with possible poles at x n of order < d p and they still define 
an integrable system on V a . In particular, h P2 takes values in H°(K 2 (2 x n )) which 
is dual to H 1 (K~ 1 (—2J2x n )), the tangent space to the moduli space of curves E with 
fixed punctures x n and first jets at x n of holomorphic local parameters z n , z n (x n ) = 0. 
The corresponding Beltrami differentials 5/i behave like 0(z%) around x n and they are 
taken modulo d{5^) where the vector fields £ are also O(z^) around x n (such vector fields 
do not change the first jets at x n of the local parameters z n ). 5fi may still be coupled 
to quadratic differentials G H°(K 2 (2 £ x n )) by (JlT^) and Eq. (|L9|) defines for [6/j] G 
H 1 (K^ 1 (-2J2 x n)) Hamiltonians on V that are in involution. 



2 Knizhnik-Zamolodchikov-Bernard connection 

The phase space V = T*M may be (geometrically) quantized by considering the space 
H°(C k ) of holomorphic sections of the k th power of the determinant line bundle £ over 
J\f (more exactly, over its semi-stable version Af ss ) as the space of quantum states. Such 
sections are given by holomorphic functions ip on A satisfying the Ward identity 

iP(A) =e- kS{h ' A ^{ h ~ 1 A) (2.1) 

for h G Q and with S(h, A) standing for the action of the gauged Wess-Zumino-Novikov- 
Witten (WZNW) model. The identity (^]l]) expresses the gauge invariance on the quantum 
level. The vector spaces H°(C k ) arise naturally in the context of the WZNW model and 



of the Chern-Simons theory [25]. They are finite-dimensional and their dimension is given 
by the Verlinde formula [24]. Put together for different complex structures of E, they form 
a holomorphic vector bundle W over the moduli space of complex curves. In the language 
of functions ip, the <9-operator of this bundle is given by 

M = ( d ¥ + Wi / E tr A) * (2-2) 
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where dj^ differentiates ip viewed as a function of the unitary gauge field B = —A* + A = 
—A\dz + A s dz (functions of B are naturally identified for different complex structures 
on £). The bundle W may be equipped with a projectively flat connection V KZB [25] Jl7| 



which may be traced back to the works of Knizhnik-Zamolodchikov [18] and Bernard 



In the present description of W, the KZB connection takes the form [12] 

vJ^V = (*„ - J^iix&ti ~ t J^:m(m^)l)iP, (2-3) 

V| ZB ^ = %V- (2-4) 

where n = k+g v with g v denoting the dual Coxeter number of G. The symbol * * indicates 
that one should remove the singularity at the coinciding points of sl(y) ^ before setting 
x = y. How this is precisely done depends on some choices (e.g. of a projective connection 
or a metric on each E) but the choices lead to connections differing by addition of a scalar 
form. 



The second order operator on the right hand side of Eq. ( |2.3| ) has the principal symbol 
(obtained by replacement of Jj by proportional to the Hitchin Hamiltonian hg^. 

The KZB connection V^ ZB may be considered a quantization of ^hs^ which, instead of 
acting in a fixed space H°(C k ) relates two such spaces for the complex structures differing 



by <fyi[17|. Note that if we rescale 5fi i— > nSfi, we should obtain from the KZB connection in 
the limit k — > operators acting in the space H°(C~ 9 ) corresponding to a fixed complex 
structure. This space becomes non-trivial if we admit singular sections of C or work with 
higher cohomologies of C~ 9 . It also admits a quantization of the non-quadratic Hitchin 
Hamiltonians [||[11]. For k ^ — g w we may also obtain from Eq. (2.3) operators in a single 



space if we chose a local trivialization of the bundle W (or of a bundle W D W). 

The above quantization extends to the case of the phase space V if the coadjoint orbits 
O n associated to points x n G £ correspond to irreducible holomorphic representations of 
G in vector spaces V n (i.e. to irreducible unitary representations of the compact form of 
G). The quantum states are now represented by holomorphic maps on A with values in 
V = <g> V n satisfying the Ward identities 

n 

ip(A) = e~ k 5 ^' A ) <g> h(x n ) ip{ h ~A) (2.5) 

n 

for h G Q generalizing Eq. ( |2.1[ ). The spaces of solutions are still finite-dimensional and 
form a holomorphic vector bundle over the moduli space of punctured curves with first jets 
of local parameters at the punctures. The complex structure and the KZB connection are 
given by the same formulae with Beltrami differentials 5(j> restricted to behave like O(z^) 

jZ^ip = 0{z~ l ) and tr ; & £ {x) & ^ \ 

Sp, 

, , , onrl oil f no ohn 



at the punctures. Since jjbx ip = 0{z^ 1 ) and ^l sA(x) 5A(x) » ^ = ®( z n 2 ) around x n , there 
is no problem of convergence of the integrals over S. Again, V^p B may be viewed as the 
quantization of the Hitchin Hamiltonian Tp-^hfin and all the above remarks apply. 
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3 Genus zero 



Up to diffeomorphisms, there is only one Riemann surface of genus zero: the Riemann 
sphere P 1 = Cu {oo}. On P , the gauge orbit of the zero gauge field is open and dense in 
A, i.e. the generic gauge field takes the form 

A = h~ l dh (3.1) 

where h E Q is determined up to left multiplication by a constant g £ G. The equation 
n(A, <£) = J2n^n8 Zn , with X n belonging to the (co)adjoint orbit O n associated to the 
puncture z n , becomes 

d( h $) = J2"n6 Zn (3.2) 

n 

where v n = h(z n )\ n h(z n )~ l € O n . This equation has a (unique) solution 

^) = E^r£ (3-3) 



if and only if the sum of residues is zero, i.e. if Yin v n = 0- We obtain then for V D defined 
by Eq. ( flTi] ): 

V = {uexO n \j2^n = 0}/G. (3.4) 

The (co) adjoint orbits carry a natural symplectic structure leading to the Poisson bracket 
{v a ,v b } = if abc v b for v a = trt a v where t a are the generators of Lie(G) s.t. tr t a t b = 5 ab 
and [i a ,i 6 ] = if abc t c . It is easy to check that the complex symplectic structure on V 
coincides with the one obtained by the symplectic reduction of x n O n with respect to the 
diagonal action of G. 

The Hamiltonians h p are 

n Z Zn 

In the special case p2 = ^tr, 
where 



<5„ = i«, h n =J2^^- (3-7) 



Zi, 



Let 5fj, be a Beltrami differential regular at infinity and behaving like 0((z — z n ) 2 ) around 
the insertions. Necessarily, 5fi = <9(<5£) for a regular vector field 5£, = 5^ z d z on P . <5£ 
is determined up to infinitesimal Mobius transformations (a + bz + cz 2 )d z . We may take 
z 1 = z + <5£ z as the new complex coordinate on P 1 with the modified complex structure. 
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The modification is then equivalent to the shift 5z n = 5£ z (z n ) of the insertion points 
and the shift b\n = Xnd z (5^)(z n ) of the first jet of the local parameter at the punctures 
parametrized by the 9 2 -derivative Xn of the parameter at z n . An easy calculation involving 
cutting out small balls around the insertions and integration by parts shows that 

h *K = E/ ( ( 1 S2 5 n + ^~ O(g-) 2 ^) 
^ JY,\{Z - Z n ) 1 Z-Z n J 2m 

= ~hi E ( 6n Xn l6 Xn + K 5z n ) . (3.8) 



The quantum states if) at genus zero may be labeled by their values V'(O) at ^4 = 
which belong to the subspace V G of the G-invariant tensors in V = <8> n Vn- Indeed, ip is 
determined by its values on the dense (/-orbit of A = 0. Hence the bundle W is a subbundle 
of the trivial bundle with the fiber V G . The KZ(B) connection reduces in this case [^] to 
the formula 

VgV(O) = ( 6 Xn(d Xn - Xn'An) + Sz n (d Zn - l -H n )\ m , (3.9) 

V f^(°) = E (^An + ^(0) (3.10) 

n 

where 

An=2K"in*n; #n = E ( 3 -H) 

with denoting the action of the generator t a in the factor V n of V. A n is a number, the 
conformal weight assigned in the WZNW theory to the irreducible representation of G in 
V n Note that, modulo the shift i— > k = k + g v , A n and ^-Hn may be obtained from 
(2ni)' 2 ^ n an< ^ finTp res P ec ti ve ly) by the (geometric) quantization of the coadjoint orbits 
which replaces the functions by the operators so that the Poisson bracket turns 
into M times the commutator (% plays the role of the Planck constant). The flatness of 
the connection V KZ , (V KZ ) 2 = 0, follows from the equation [H n ,H m ] = 0, equivalent to 

t a t a 

the classical Yang-Baxter equation (CYBE) for - L - 2 -: 



j.a j.a fO, fa 



+ cyclic permutations = 0. (3-12) 



4 Genus one 

At genus one, every Riemann surface is isomorphic to an elliptic curve E T = C/(Z + rZ) 
where r is a complex number of positive imaginary part ti. Denote by A (A+) the set of 
(positive) roots of Lie(G), by e a the step generators attached to the roots a and let (rf) 
be an orthonormal basis of the Cartan algebra Lie(T). We set u a = trua and u 3 = trurf 
for u E Lie(T). We shall need some elliptic functions: the Jacobi theta function 

oo 
£=-oo 
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the Green function P x of the twisted <9-operator 

with the properties P x (z + 1) = P E (z), + r) = e" 2 ™^^) and P x (z) = \ + 0(1) 
around z = 0, the function 

P = #iM (4.3) 

s.t. p(z + 1) = p(z), p{z + r) = p(z) — 2-iri, p{z) = \ + 0(1) around z = and, finally, the 
Weyl-Kac denominator 

oo 

n(u) = e 2 ^ Td / 24 JJ (e** 40 - e"™ Q ) J | [(1 - e 27ri£r ) r [] (1 - e 2 ™ £r e 2 ™")] (4.4) 

aeA + 1=1 oeA 

where d denotes the dimension and r the rank of G. 

On E T a generic gauge field is in the orbit of A u = irudz/T2, for u € Lie(T), i.e. 

A = ^k" 1 = (h u h)~ l B(h u h) (4.5) 

where /i u = e n ( uz - uz )/ T 2 . Consequently, the gauge fields may be parametrized by u and h. 
To avoid ambiguities, we have to identify the pairs as follows 

(u,h) ~ (wuw' 1 , wh) ~ (u + g v , /i~v /i) ~ (u + rg v , h~^ v h), (4.6) 

for g v in the coroot lattice Q y and ^ in the the normalizer N of Lie(T) in G. Similarly to 
the genus zero case, we have to solve the equation 

QpHuh^ = Y,v n 5 Zn (4.7) 

n 

where v n = (h u h)(z n ) X n (h u h)(z n )~ l ■ Decomposing v n = J2 a u n ae a + ^n witn v n = v hf £ 
Lie(T), we can solve the above equation if and only if J2n v n = 0- I n that case, 



t (4-8) 



huh Hz) = Lo + T,(E P uM-Zn)^ a e a + p(z-z n )v°)\ 

\ n a / 

for an arbitrary constant ifQ = (f^rf . Performing the symplectic reduction, we find 

V ~ {(u,w,K) £T*Lie(T) x (xO n )\Y,v°n = 0} / Nk(Q v +tQ v ) (4.9) 

n ' 

where the action of Nx\(Q v + tQ w ) implements the identifications 

(u, (fo,K) ~ (wuw' 1 , wipow' 1 , wvw~ l ) ~ (u + q y , y?o, {h~v{z n )v n h q v(z n 
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The symplectic structure of V is that of the reduction of T*Lie(T) x (x O n ) by the group 

n 

iVxi(Q v + rQ v ). Now it is easy to write down the Hitchin Hamiltonians. Let us us do it 
for p2 = ^tr. A straightforward computation identifying the pole terms leads to 

h P2 = {~^P\ Z ~ z n)$n + ^p{z - z n )h n + h^ij^j (4.10) 

n n 

where, as before, 5 n = hv%u® an d 



2"n"n 

r r a// 

h Q = \H + \ E { E d x P u a {z n - Z rn )uy- a + | E ~jt( Zn - ^K^m}' O 4 - 11 ) 

j=l m,n a j=l 1 



K = E V n^l + E ( E ( z n ~ Zm) V> m ° + E P^ ~ z m) ■ ( 4 - 12 ) 

3=1 m=£n a j=l 



Note the similarity to the genus case (3.6), 

Let 5fi = 8\i%d z dz be a Beltrami differential on E T behaving like 0((z — z n ) 2 ) around 
the insertions. The modified complex structure corresponds to the complex coordinate 
z > = z + §^£t + 5i z s.t. d z z' = Sfi%. We require that 5£ z (z + 1) = 8£ z {z + r) = 8£ z (z). 
5t is determined from the condition that the integral of <5/i§ over E T is equal to that of 
2^-5r. 5^ z is unique up to an additive constant. Note that z'(z + 1) = z'(z) + 1 whereas 
z {z + r) = z' + r' where t' = r + Jr. Hence the deformed curve is isomorphic to E T > with 
the punctures moved to z' n = z n + 5z n and the first jets of local parameters changed to 

Xn = Xn + SXn with 



5 Zn = ^-^St + 8e{z n ) , Xn^Xn = ^ + 9M z {z n ) . (4.13) 



Again by a straightforward calculation substituting 8 fig = dgz', cutting out small balls 
around points z n and integrating by parts, we obtain 



h P2 5fi = ^E (^Xn^Xn + h n 8z n + ^-h 5T^j . (4.14) 



The quantum states ip at genus one may be characterized by giving holomorphic func- 
tions ip{u) on Lie(T) with values in V T , the subspace of T- invariant tensors in the product 
V of the representation spaces, 

i>.{u) = n(u)e-^ tru2/(2T2) ® ( e -^n-z n )ulr 2 \ ^ Au y ( 415 ) 

n V In 

The KZB connection takes the form [||, [5], [l0|, [|, || 

V& ZB ^ = E fcn(9 x „ - X^A n ) + 8z n (d Zn - l -H n ) + 5r(d T - ^H )) (4.16) 

vf 5 ^ = E (^An + ^9*. + E^) <a (4.17) 

n 



8 



where A n is as before and the operators \Hq and \H n are obtained from the Hamiltonians 



Tjho and (9 h n by the replacement 



(2m)' 2 u (27ri) 



'o ' * S"e Qn , y? n i-» ^rf n , (4.18) 



i.e. by the geometric quantization. The resulting Hamiltonians Hq and H n act on general 
meromorphic functions on Lie(T) with values in V T . The flatness of the KZB connection 
is ensured by their commutation: [H n , H m ] = 0, for n, m = 0, 1, . . ., following from the so 
called dynamical CYBE M. 



5 Genus two 



5.1 Curve and its Jacobian 

Let S be a curve of genus 2. Choosing a marking, i.e. a symplectic homology basis 
(A a ,B a ), a = 1,2, on S, we may fix the corresponding basis (co a ) of the holomorphic 
1,0-forms (abelian differentials) s.t. J^ a uj b = 5 ab . The i? a -periods of the abelian differen- 
tials give rise to the symmetric period matrix r, r ab = J Ba uj b , with a positive imaginary 
part. The map 

z{x) = (5-1) 

realizes S as a double covering of P 1 ramified over six Weierstrass points x n or as a 
hyperelliptic curve given by the equation 

6 

y 2 =Y[(z-z n ). (5.2) 

n=l 

The coordinates z n = z(x n ) are assumed to be finite. This may be always achieved by 
an appropriate choice of the marking. Curve X may be viewed as composed of the points 
(y, z) and of two points at infinity. The covering of P 1 is (y, z) i— > z. In this representation, 
the holomorphic 1,0-forms and the holomorphic quadratic differentials on S have the form 

to = (a + bz)dz/y , (3 = (a + bz + cz 2 )(dz) 2 /y 2 , (5.3) 



respectively. In particular, 

uj 1 oc dz/y, uJ 2 oc zdz/y (5-4) 

with the same proportionality constant. 

The Jacobian J 1 of the degree 1 holomorphic line bundles may be represented as 
C 2 /(Z 2 + rZ 2 ). In particular, the spin structures L, L 2 = K, correspond to points e + re' 
with e G iZ 2 /Z 2 . There are 6 odd spin structures e n + re' n labeled by the Weierstrass 
points z n , the zeros of the holomorphic sections they admit, and 10 even spin structures 
without holomorphic sections. 
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5.2 Theta functions 

The degree 1 bundles with holomorphic sections form the theta-divisor in J 1 . The holomor- 
phic sections of the fc th -power of the corresponding theta bundle over J 1 may be represented 
by the holomorphic theta functions of degree k on C 2 defined by the relations 

e™ kn ' Tn+2 ™ kn - u 9(u + m + Tn) = 6(u) (5.5) 

for m,n G Z 2 . They form the space ©& of dimension k 2 . For k = 1 there is a single (up to 
normalization) theta function 

0(u) = J2 e™ n - Tn+2 ™ n - u , (5.6) 

rial? 

the Riemann theta function. For k = 2 there are four independent ^-functions. One can 
take them as 

Q e (u} = e 2ni ( n + e )- r ( n + e ) + 47ri (n+e)-u (5 7) 



res 



for e G ±Z 2 /Z 2 . 

i?(u + - u) = ^ e (u)# e (v) (5.8) 

e 

is a second order theta function in both u and i>. The map v i— ► — u) determines 

an embedding of the Kummer surface J 1 /^ = C 2 /(Z 2 + rZ 2 )/Z2 onto a quartic surface 
/C in the 3-dimensional projective space P02 (Z2 maps the degree 1 line bundles L into 
L~ 1 K or v eC 2 into — v). 



The double theta function Q5.8] ) determines a non-degenerate symmetric quadratic form 
on the space 02 dual to 02- It permits to identify ©2 with ©2 by sending <j> G ©2 to 
l{4>) G ©2 defined by 

i(<t>)(u) = + ■)>#>• ( 5 - 9 ) 

The identification exchanges the basis (9 e ) of ©2 with the dual basis (0*) and the Kummer 
quartic /C with its dual version JC* C P@2- JC* is composed of linear forms proportional to 
the evaluation forms (j) u defined by 

(9,<f> u } = 9(u). (5.10) 

The group (^-Z/Z) 4 of spin structures acts on ©& for even k by endomorphisms U ee > 
defined by 

(U e ^6)(u) = e wike '- Te ' +2wike '- u 6(u + e + Te'). (5.11) 

For k not divisible by 4 this action is only projective: U ee >Uf f> = (— l) 4e '-^ U e+ f j( ,' + f and 
it lifts to a Heisenberg group. For k = 2, 

U e:e '6 e ii = ( — l) 4e ' e 6 e i +e ii . (5.12) 

The action of U ee i preserves the Kummer quartic JC C P@2 and the action of the transposed 
endomorphisms U J , preserves JC*. 
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5.3 Moduli space of SX 2 -bundles 



In the fundamental paper [21], Narasimhan and Ramanan proved that the moduli space 
M s of the stable SL2 holomorphic bundles is canonically isomorphic to PO2 \ JC. The 
isomorphism associates to an SX 2 -bundle E the second order theta function 9 vanishing at 
the points u G C 2 corresponding to the duals of the line-subbundles of the rank 2 bundle 
associated to E. In other words, if A is the gauge field whose £?-orbit corresponds to E 
and if L u denotes the degree 1 line bundle corresponding to u G C 2 /(Z 2 + rZ 2 ) then the 
theta function 9 associated to E vanishes at u if and only if there exists a pair s = (si, s 2 ) 
composed of sections of L u s.t. (d + A)s = 0. The semistable compactification of the 
moduli space M s is 

Mss = Pe 2 (5.13) 

and, exceptionally, it is smooth. The points of the Kummer quartic /C represent (classes 
of) the semistable but not stable bundles. Hence for G = SL2 the phase space of the 
Hitchin system on the genus 2 curve with no insertions is 

t*m ss t*pb 2 ^ { (e, <f>) g e 2 x e 2 * \e^o, (#,</>) = o}/c x (5.14) 

with the action of t G C x given by (0,4>) i— » (t9, t^ 1 ^). As a symplectic space, it is the 
symplectic reduction of T*(0 2 \ {0}) by the action of C x . Using the bases (9 e ) and (9*) 
to decompose 

9 = Y / QeO e , <f> = J2Pe0*e, (5-15) 

we may represent T*M SS as the space of pairs (q,p) G C 4 x C 4 , q ^ 0, q ■ p = 0, with the 
identification (q,p) ~ (tq, t~ 1 p) and the symplectic form induced from dp ■ dq. 

5.4 Hitchin map for G = SL 2 

The Hitchin map h P2 : T*P0 2 — > H°(K 2 ) appears to take a particularly simple form 
resembling the genus formula (|3.6| ): 

^ = ;f - — ^ — -(^) 2 = e-^(^) 2 (s-ie) 

P2 2 n ^ (z-z n )(z-z m ) v 2-^ ^ z _ Zn ^ m > v ; 



where 



r nm (9, <t>) = ± (U en< 9, U e l e ,J) (U emKi 9, U^J) (5.17) 
(the last two factors on the right hand side coincide modulo sign) and 

K = £ . (5.18) 



With the help of Eq. (|5.12|) , r nm 's may be rewritten in the language of g's and p's as explicit 



homogeneous polynomials of order 2 in q e and in p e , see below. The identity 

E r — = ( 5 - 19 ) 
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holding for each n guarantees that T\(z — z n ) Y] , _ r \P l _ — r is a quadratic polynomial 

n =l n,m=l (Z Zn>[Z Zm> 

in z. It follows that the right hand side of Eq. ( 5.16j ) determines a quadratic differential 



on £ of the general form given by Eq. ( [5.31 ). 

The equality ( |5. 16 ) is not immediate. It was established in four steps. We shall only 



enumerate them here. The two first crucial steps were performed in [14|. It was shown 
there that for any 9^0 and u G C 2 s.t. 9{u) = 0, 

MMJ = -^{duaeiu)^) 2 . (5.20) 

The above equation describes the quadratic polynomial h P2 (9, ■) (with values in H°(K 2 )) 
on the quartic K* e = JC* n in the projectivized subspace of ©2 perpendicular to 9. 
In principal, it determines h P2 completely. It was observed then that the above formula 
implies that for each Weierstrass point z n the conic 

C n = {C x G F9 1 | MM)k = 0} (5.21) 

is in fact the union of two bitangents to JCq . The explicit equations for the bitangents to the 
Kummer quartics known since about a century permitted then to establish Eq. ( |5.16j ) up 



to multiplication by a ^-dependent factor [14|. The other steps in the proof of the formula 



( 5.16| ) were taken in |l3[] were it was established that the Hitchin map h P2 possesses the 



important self-duality property: 

h P2 r l (6)) = MM) (5-22) 

or that h P2 (q,p) = h p2 (p,q) in the language of ( ^.15| ). This property, far from obvious in 
the original formulation of the Hitchin system, restricted the ambiguity on the right hand 



side of Eq. (5.16) to a (possibly curve-dependent) constant factor. The latter was fixed in 
[|l3| by a tedious calculation of h P2 at special points of K, x JC* . 



5.5 Deformations of complex structure 

For the hyperelliptic curve, the variations of the complex structure described by the Bel- 
trami differentials on S change the images z n of the ramification points of the covering 
of P 1 . Let us find these changes. Let 

J a = uj a + 5oj a + 5uj a (5.23) 

denote a deformed basis (u/ a ) of the abelian differentials with 5u) a of 1,0-type and 5uj a of 
0,1-type in the original complex structure. 5u) a and 5u) a have to satisfy the relations 

5uj a = oj a 5fi, 3{5uj a ) = -d{oj a 8fi) (5.24) 

stating, respectively, that u>' a is of the 1,0-type in the deformed structure and that it 
is a closed form. The equation for Sco a always has solutions. They are defined modulo 
abelian differentials and the ambiguity may be fixed by demanding that J Aa u>' b = 5 ab . The 
deformed covering map onto P 1 is then 

II \ Oj' 2 (x) 50J 2 x^ 1 / \ lr 
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The ramification points x' n of the map z' are determined by solving the equation 

&z'{x' n ) = (5.26) 

which, upon rewriting x' n = x n + 5x n becomes 

x, .2 r. .1 

(d) 2 z(x n )5x n + d(— ){x n ) -d(z— ){x n ) = 0. (5.27) 

Since dz(x n ) = 0, the quadratic differential (d) 2 z(x n ) is well defined. Besides, since 
the ramification points are isolated, it does not vanish so that one may solve the above 
equations for 5x n . The ramification points x' n are mapped to 

5uj 2 8oj^ 

z' n = z(x n + 5x n ) = z n H 7-{x n ) - z n — r-(x n ). (5.28) 

uj l to 1 

We infer that 

8co 2 Su)^ 
5z n = z' n - z n = — r (x n ) - z n —-(x n ) (5.29) 

are the variations of z n corresponding to the Beltrami differential e)/i. 

We may now find the values of the Hitchin Hamiltonians h$ u = J s h P2 5fi related to 
the Beltrami differentials. Note that, by virtue of Eqs. ( |5.25| ) and ( |5.24| ), 

d(8z) = d{z' -z) = diSuo 2 )/^ 1 - zBidio 1 )/^ 1 = -d^S^/uo 1 

+ zd(uj 1 5fi)/uj 1 = -diz^d^/to 1 + zd^S^/cu 1 = (dz)dfi. (5.30) 

A straightforward integration by parts over the region in £ without small balls around the 
Weierstrass points x n and around 2 points at infinity gives now: 



\-^K 5z n . (5.31) 



7T2 

n=l 



The comparison with Eq. ( |3.8[) shows an additional factor 2 which comes from the double 
covering. 

5.6 Relation to Neumann systems 

The Poisson-commutation of the Hamiltonians h n (of which any 3 give independent action 
variables of our integrable system) is equivalent to the CYBdike equations 

{ ^ , -J^} + cycl. = 0, | rnrn ^ ^J^ } = Q (g 3^ 

for {re, m} n {p, q} = 0. The above relations may may be directly checked, as noticed 



in ||i4|| . The more recent observations of the paper |15] on the Knizhnik-Zamolodchikov 



Bernard connection in the same setup (see below) permit to identify the integrable system 
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with the Hamiltonians h n of Eq. ( |5.18| ): it is a modified version of the classical genus 2 
Neumann systems [^, |l^, || whose original version is also rooted in the modular geometry 
of hyperelliptic curves. This goes as follows. 

The phase space T*P 3 (without the zero section) may be identified with the coadjoint 
orbit (D\ of the complex group SX4 composed of the traceless rank 1 matrices The 
action of SL4 in A 2 C 4 preserves the quadratic form induced by the exterior product on 
A 2 C 4 and the identification A 4 C 4 = C 1 . It leads to the double covering SX4 — > SOq if we 
choose in A 2 C 4 the Pliicker basis turning the quadratic form into the sum of squares. The 
inverse relation is the complexified version of the twistor calculus. Upon the identification 
of sU with so6, the coadjoint orbit 0\ becomes the one composed of (complex) rank 2 
antisymmetric matrices J = (J n m) of square zero: J 2 = 0. Such matrices are of the form 
Jnm = QnPm — PnQm with vectors P, Q € C 6 spanning an isotropic subspace, i.e. with 
Q 2 = Q ■ P = P 2 = 0. Given (q,p) with q ■ p = 0, in order to find the corresponding pair 
(Q, P), it is enough to complete vector q to a basis (q, f±, f 2 , fy) of C 4 s.t. f\ -p = f 2 -p = 
and /3 • p = 1 and to set 



Q = qAf U R = f 2 Af 3 , P = qAf 2 /Q-R 
in the language of A 2 C 4 = C 6 . In the Pliicker coordinates, we may take 

Qi = -(qipz + qm) , Q2 = i(qiP3 - qm) , Q3 = -i(qip 2 + qm) , 



(5.33) 



Q4 
Pi 



q\P2 - qm , 

_1 g2P4+<?3Pl 
2 <?2P2+<?3P3 ' 

1 <?3P4~g2Pl 

2 <J2P2+I?3P3 ' 



Q 5 
P 2 
P 5 



q2P2 + qm , 



l_ <?2P4~g3Pl 
2 (J2P2+93P3 

1 

2 ' 



Q 6 
Pz 
P(i 



i(q2P2 + <?3P3) , 



i_ <?3P4+g2Pi 

2 <?2P2+<?3P3 



where q x = q {0fi) , q 2 



m0)' 



qi 



9(0 i)> 94 — 3(1 I) an d similarly for p's. In terms of 



(Q, P), the symplectic form is dP ■ dQ. The functions J nm on 0\ have the Poisson bracket 



{Jnmi Jmp} 
{Jnmi Jpq\ 



J, 



np 







A straightforward check shows now that 
J12 = |(<?iPl + q2P2 ~ q3P3 ~ qAPA) , 

^(qiPA + <?2P3 - 93P2 - qAPl) , 
^(<?l£>3 + <?2P4 + q3Pl + 94P2) , 
~\{q\P4 + <M?3 + 93P2 + q^pi ) , 



Jl4 
.^16 
«/24 

<^26 
^35 
<-?45 
<^56 



9 (<7lP3 + 92P4 - q3Pl ~ qiP2) 



\{q\P2 + q%P\ + 93P4 + 94P3) 
\{q\P2 - q2P\ + 93P4 - 94P3) 
f (<M>1 - 92P2 - q3P3 + 94P4) 



for n,m,p different, (5.34) 
for n,m,p,q different (5.35) 



J13 = -^(qiP4:-q2P3-q3P2-qm), 

J15 = -\{qiP3 - q2P± - q3Pi + q±P2) , 

J23 = -\{qip± - q2P3 + q3P2 - qm) , 

J25 = ^(qiP3 - q2P4 + q3Pi - q^P2) , (5.36) 

J34 = ~^{qiPi - q2P2 + q3P3 - <M>4) , 

J36 = ~^(qiP2 - q2P\ - q3Pi + q±P3) , 

J 4 e = -\{q\P2 + q2P\ - q3Pi - qAP3) , 
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and that 



Trim — a (Jwn) • (5.37) 

The equations (15.321) assuring that the Hamiltonians h n = — -? V Poisson-commute 

^ * z n z m 

m^tn 



follow directly from the soq algebra ( 5.34| ), ( |5.35| ). The original Neumann systems are very 



similar but involve the coadjoint orbits of SOn composed from rank 2 antisymmetric 
matrices of square ^ ||. Such orbits, contrary to the one we consider, have nontrivial 
standard real forms. 

5.7 Lax matrix approach 

Although the change of the orbit modifies dimensional counts and many details, the meth- 
ods used in the analysis of the Neumann systems, in particular the Lax method developed 
in jjj, generalize with minor variations to our system and permit to find explicitly the angle 
variables of the genus 2 Hitchin system. The Lax matrix may be taken as L(() = (L nm (£)) 
with 

LnrniC) = C Jnm Z n 5 nm . (5.38) 

As in [0, the Poisson brackets (|5.34Q , (|5.35|) may be rewritten in the matrix form as 

{L(C)®1,1®L(C0} = [£(C)®l,r-(C,C0] - [l®L(C'),r + (CO] (5.39) 
where the r-matrices 

r^CCO = ^rC ±^T (5.40) 

with C mn . 

qp — dmqfinp and T m n,qp — ^mp^nq satisfy the CYBE. The above form of the 
Poisson bracket implies immediately that 



{trL(C/, trL(C') £ '} = (5.41) 



for all C) C- Since 



I J 2 

the Hamiltonians h n = — j - may be expressed as combinations of the quantities 
trL((^Y. It is not difficult to see that the converse is also true. 
More generally, Eq. ( [5.40 ) implies that 



{trL(C) £ ,L(C')} = [M e (C,C'),L(C')] (5.43) 

with 

M e (C, CO = igr + i ^ L (-0 £_1 • (5-44) 
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It follows that the commuting time evolutions of the Lax matrix L(C) generated by the 
Hamiltonians h n , 

S n L(0 = {h n ,L(Q}6t n , (5.45) 

are isospectral. In other words, the spectral curve S given by the characteristic equation 

det(L(C) -z) = (5.46) 

is left invariant by the dynamics generated by any of the Hamiltonians h n . An easy 
calculation using the fact that the matrix J has rank 2 gives 



det(L(C)-,) = J] (*-*»)(! + fC 2 £ U ft ~ - ( 5 - 47 ) 
Upon the substitution a = ^ Y\ n (z — z n ), the characteristic equation ( |5.46 ) becomes 

o 2 = \ mz-z n f e (z _ z ir z - z ) 55 (5 - 48) 

Since P(z) is a polynomial in z of order 8 (see the remark after Eq. ( |5.19| )), this is the 
equation of a hyperelliptic curve S of genus 3 composed of pairs (a, z) and 2 points 
corresponding to z = oo. We shall consider only such points of the phase space that S is 
smooth. The 1,0-forms 

n b = z b dz/a (5.49) 
with b = 0, 1, 2 form a basis of the abelian differentials on S. 

We shall search for the eigenvectors X = (X n ) of the Lax matrix. This will allow to 
adapt the arguments described in great detail in Sect. 4 of [19] to the present case. The 
eigenvector equations 

CJ nm X m = (Q n (P-X)-(P n (Q-X) = (z-z n )X n (5.50) 

imply that 

[z-z n )X n = aQ n + bP n . (5.51) 

Upon multiplication by a = rC^Iln^ — z n)i Eq. ( [5.50 ) becomes a system of 2 linear 
equations for a and b: 

(a + V)a + iWb = 0, 

-iUa + (a-V)b = (5.52) 

where 

Qn 



u(z) = n 



n=l n=l " 
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6 6 n p 

v(z) = *!!(*-*»)£ 



n=l n=l Z Zn 

6 6 p2 



W(z) = l[( z - Zn )Y^^^ (5.53) 

n=l n=l n 

are 4 th -order polynomials in z. The non-trivial solution exists if 

a 2 = U{z)W{z) + V(z) 2 = P{z) (5.54) 

where the last equality follows by a straightforward check. The system ( 5.52| ) of linear 
equations defines a holomorphic line subbundle L of the rank 2 bundle W = C 2 (g> 0(4p+ + 
^Poo) over t ne hyperelliptic curve 5 (the coefficients behave as z 4 at infinity). As solutions 
of (|5,52j ) we may take for example 

a = a-V(z), b = iU(z) or a = -iW(z) , b = a + V{z). (5.55) 

Since a and b are proportional to z 4 at infinity, they define holomorphic sections of L C W 
regular at p^. They vanish at four points 

p' a = (V(z' a ),z' a ) or P'L = (-V(4),4), (5.56) 

respectively, where z' a are the roots of U and z" a are those of W. Hence the degree of the 
line bundle L is equal to 4. H°(L) has dimension 2 and is spanned by the two solutions 

(HI). 

5.8 Angle variables 

The knowledge of the bundle L may be encoded in the image w of L in the Jacobian J 4 (5) 
under the Abel map: 

4 Pa 4 Pa 

W = J2jn = J2jn (5.57) 

a=1 Po a=1 po 

were fi = (f2 b ) is the vector of the abelian differentials ( 5.49| ) on S and po is a fixed point 
of S. Under the infinitesimal time evolution ( |5.45| ) inducing the changes 5 n U ', <5 n y, 5 n W, 
of the polynomials U,V,W, the image of the Abel map changes by 

4 / b r / 4 lib j- a 

^ -£><^- (5 ' 58> 

The variations of the zeros of f7 are 

^ - " iT^y (5 - 59) 

and similarly for 8 n z'^. A direct calculation gives 
S n U(z) = {h n , U(z)}5t n = 4< J] (*» ~ a™)" 1 £^1^ ^ , (5.60) 
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see |l9t| , page 3.69. Hence 



4 / b 

5 n w b = 4iY[(z n - zm)- 1 U(z n ) 5t n ]T ( , _\ jv( , \ ■ ( 5 - 61 ) 

m+n a=l ^ Z ™> U \ Z a) 

The vanishing of the sum of residues of the meromorphic form j~r z vf/ra hnplies that the 
last sum is equal to — z h n U{z n )~ 1 so that 

5w b = {w b ,h n }5t n = j l[ (z n - z m )- 1 z b n 5t n (5.62) 

which does not depend on the phase-space variables. We infer that the Hamiltonians h n 
generate constant flows on the complex torus J 4 (5) = C 3 /A where A is the lattice of 
periods of Q. Modulo a constant linear transformation, the coordinates w b , b = 0,1,2, 
provide together with 3 of the Hamiltonians h n a Darboux coordinate system for T*P 3 . 
We have thus found the angle variables of the Hitchin system (as the angles of J 4 (5)). 

It should be stressed that the above approach based on the Lax matrix is simpler then 
the one obtained by following the general procedure for the Hitchin systems. In particular, 
Eq. ( |1.7| ) giving the spectral curve C in the general approach is, as shown in Ji~3| ], 

e = -4nc*-*o e 

n=l n,m=l 

to which one has to add the equation (|5.2j) of the original curve S of genus 2. C is of genus 5 
and it is a ramified cover of both S (by forgetting £) and S (by setting a = j^Cl/)- While 
the general construction would give the angle variables as those of a 3-dimensional Prym 
variety in the 5-dimensional Jacobian of degree —2 line bundles on C, the Lax approach 
gave them as the angles of the degree 4 Jacobian of S. 



% _ v (5-63) 



5.9 KZB connection 

The determinant bundle C over the moduli space J\f ss = PO2 of the holomorphic SL2- 
bundles over the genus 2 curve coincides with the dual of the tautological bundle on P©2 
so that H°(C k ) is the space of homogeneous polynomials \P of degree k on the space 02- The 
Lie algebra soq acts in the space 02 by the first order differential operators, still denoted by 
J nm , satisfying the commutation relations Q5.34 ), ( 5.35| ) with the Poisson bracket replaced 



by the commutator. In the (p, g)-language they are obtained by replacing p n 's in the 
expressions ( |5.36[ ) for J nm by —d Xn . The KZB connection for the case in question has been 
work out in [jlq]. It takes the form (up to a scalar 1-form) 

Vf M ZB * = (d Zn - l -H n )V, (5.64) 

n=l 
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where 

Hn = ~\ £ (5.66) 

so that \H n is a quantization of T^fw^n obtained by the replacement J nm \— > j^Jnm in the 
classical expression for h n . The quantum Hamiltonians H n , n = 1, ... ,6, are commuting 
second order differential operators on 62 — C 4 . 

6 Conclusions 

We have described above in explicit terms the Hitchin integrable systems and the Knizhnik- 
Zamolodchikov-Bernard connection in the genus 0, 1 and 2 geometries, the last case only 
for G = SL2 and with no punctures. The main original contribution of the paper is the 
construction of the angle variables of the genus 2 system. It is a modification of a similar 
construction, based on the use of a Lax matrix, for the classical Neumann system. The 
diagonalization of the quantized Hamiltonians H n entering the genus 2 KZB connection 
for G = SL2 as well as the identification of the genus 2 Hitchin systems with punctures 
and for different groups remain open problems. 
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